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Abstract 

' The dissipative corrections to the hydrodynamic equations describing the 

evolution of energy-momentum tensor and parton densities are derived in a 
simple way using the scaling approximation for the expanding quark gluon 
' plasma at finite baryon density. This procedure has been extended to study 

I the process of chemical equilibration using a set of rate equations appropriate 

for a viscous quark gluon plasma. It is found that in the presence of dissi- 
pation, the temperature of the plasma evolves slower, whereas the quark and 
gluon fugacities evolve faster than their counterparts in the ideal case without 
O ■ viscosity. 

q; PACS number(s): 25.75.Dw,12.38Mh,12.38.Bx,24.85.+p 
<D '. 
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' I- INTRODUCTION 

A phase transition from normal nuclear matter to a quark gluon plasma (QGP) is ex- 
pected in ultrarelativistic heavy ion collisions planned to be carried at the future colliders of 
RHIC and LHC If the QGP is formed during the collision, according to the usual scenario, 
it is thought that a local thermal equilibrium will set in over a typical time scale of 1 fm/c 
i]J. The plasma would then expand hydrodynamically and cool down until it reaches a 
temperature T c when the confinement phase transition starts. However, some recent studies 
have suggested that the plasma may attain kinetic equilibrium in a very short time 
scale of r = 0.2 — 0.7 fm/c, but it may not be chemically equilibrated. Earlier theoretical 
studies [[§-§] have revealed that a chemically non-equlibrated plasma cools faster than the 
equilibrated plasma, when the plasma is assumed to be ideal i.e. non dissipative. But in 
a realistic scenario the dissipative phenomena, in principle, should play an important role 
0. It needs to be mentioned here that the dissipative effects in fluid hydrodynamics can be 
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neglected when the Reynold's number (« L/X, L is the dimension of the medium and A is 
the mean free path of the constituents) is quite large. However, in nucleus-nucleus collisions, 
as the dimension of the plasma is not large compared to the mean free path of the partons, 
dissipative effects are expected to be important during the hydrodynamic expansion of the 
plasma. The formalism of dissipative hydrodynamics is, in general, complicated since the 
form of dissipative corrections depends on the choice of the rest frame |j?|^iOH . The effects of 
dissipative phenomena in quark gluon plasma have been studied in ref M, where the results 
were reported for a chemically equilibrated baryon-free plasma. It was shown that the dis- 
sipative effect enhances the entropy production leading to a dynamical path intermediate 
between idealized isentropic and isoergic ones. The energy density and temperature of the 
plasma also decrease at slower rate with the inclusion of dissipative effects. This would lead 
to greater yields of direct probes such as photons and dileptons, which are sensitive to the 
thermal history of the reaction. However, the hydrodynamic expansion for a viscous quark 
gluon plasma undergoing chemical equilibration either at zero or at finite baryon density 
has not been investigated so far. 

In this work, we study the dissipative effects in quark gluon plasma with an emphasis 
on the results at finite baryon density. We derive a set of hydrodynamical equations for 
the evolution of the energy-momentum tensor and parton densities based on simple scaling 
laws. This set of equations has been applied to study the hydrodynamical evolution of the 
viscous plasma undergoing chemical equilibration. It is found that due to viscous heating, 
the temperature evolves slower, whereas the parton fugacities evolve faster than the ideal 
plasma (without viscosity). This is contrary to the expectation that chemical equilibration 
consumes more energy which makes the plasma cool faster. 

The present work has been organized as follows. In section II, we begin with a brief 
review of the dissipative hydrodynamics and derive a set of hydrodynamical equations using 
Landau definition of rest frame. Using this set of equations, the chemical equilibration of 
a viscous baryon-free plasma is studied in section III. Finally, in section IV, we discuss 
the results and present the conclusions. In the derivations, we have used simple scaling 
law for the dissipative corrections, for which a more formal justification has been given 
in Appendix A using Boltzmann equation for the evolution of the parton densities. The 
viscosity coefficient is shown to be proportional to e/T which is also consistent with the 
finite temperature QCD calculations as shown in Appendix B. 
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II. HYDRODYNAMICS FOR A DISSIPATIVE SYSTEM 



Energy dissipation in a moving fluid occurs due to internal friction (viscosity) and heat 
exchange between different parts (thermal conduction) of the fluid. If the mean free path of 
the constituent particles is comparable to the fluid dimension, the collective quantities like 
pressure, energy density, number density, velocity etc. may vary over a distance of mean 
free path. This will lead to decay of fluid kinetic energy as heat energy. In nucleus-nucleus 
collisions, as the dimension of the plasma is not very large compared to the mean free path 
of the partons, dissipative effects may be quite important |7]]. 

In the presence of dissipation, the energy-momentum tensor T lk and number density n % 
can be decomposed into an ideal and a dissipative part |||| 

T lk = [(e + P)u i u k - Pg lk ] + r ik , (1) 

= nv? + v\ (2) 

where e, P and n are the energy density, pressure and particle number density and u l = 
7(1, v) is four- velocity in terms of local fluid velocity v(x) and r lk and v % are the dissipative 
corrections 

The relativistic hydrodynamics is based on the local conservation laws: 

diT ik = dirt = 0. (3) 

In the present context n l represents the baryon four vector and the equation din l =0 signifies 
the conservation of baryon number. However, if there is no net particle production i.e. 
chemical equilibration is complete, the above equation is applicable to the individual partons 
as well. 

The form of the dissipative terms r tk and v % depend on the definition of what constitutes 
the local rest frame of the fluid. The four velocity u l should be defined such a way that in a 
proper frame of any given fluid element, the energy and the number density are expressible in 
terms of other thermodynamic quantities by the same formulae as when dissipative processes 
are not present. It is also necessary to specify whether u l is the velocity of energy transport 
or particle transport. Accordingly, there exists two definitions for the rest frame; one due 
to Landau and other due to Eckart. In Landau approach, u % is taken as the velocity of the 
energy transport so that energy three flux T°* vanishes in a comoving frame. The Landau 
definition requires u % T %k = UiV 1 = so that T 00 = e is the energy density and n° = n is 
the number density in the proper frame of the fluid. In the Eckart definition, u % is taken as 
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the velocity of the particle transport and the particle three current, rather than the energy 
three flux vanishes in the fluid rest frame. So in the Eckart definition of rest frame, the 
particle four vector can be written as n M = (n, 0), whereas in the Landau definition of rest 
frame n M = (n, z7). Therefore, the two frames are related by a Lorentz transformation with 
a boost velocity vjn. 

It is found that due to ill defined boost velocity []7|, the energy three flux in the Eckart 
frame, which involves heat conductivity k is not well defined as it diverges in the limit 
of chemical potential fi — ►(). On the otherhand, in the Landau definition heat conduction 
enters as a correction to baryon flux. It was shown that inspite of the divergence of k, the 
correction to the baryon flux v % is finite J7| . We use the Landau definition for the subsequent 
study of the evolution of the energy momentum tensor and parton densities of the quark 
gluon plasma. 

Under the assumption that the dissipative terms r lk and v* are of first order in the 
gradients, the requirement that the entropy increases with time ( d a a a > 0, a a = au a — 
fiv a jT with chemical potential fi and temperature T ) leads to (in the Landau definition) 

m 

T a P = r]{V a u p + V p u a - - A Q/3 VX) + £A a/3 V p w p , (4) 



nT 
.e + P. 



E 

T. 



(5) 



where = g* — u a up and V a = A af3 df3, r] and £ are the co-efficients of the shear and 
bulk viscosities and k is the co-efficient of heat conduction. As mentioned before, the heat 
conductivity does not enter in the energy flux T°\ but rather as a finite baryon current 
v % oc d l (fi/T) in the rest frame of the fluid. In the presence of dissipation, the conservation 
equations (3) need to be solved with the above correction terms. Using one-dimensional 
boost-invariant scaling law |l| and Eq. (4) for r lk the identity UkdiT lh =0, resulting from the 
energy momentum conservation, takes the simple form |7| 

P 4r? + £ 



de e 
dr t 



(6) 



The evaluation of dirt = is not straightforward as one needs to know the space time depen- 
dence of /i/T and also care has to be taken so that u l is finite although k diverges in the limit 
fj, — > 0. However, we will follow a simple procedure to obtain the evolution of the baryon 
density directly without invoking the baryon current (hence the conductivity) as given by 
Eq. (5). It may be mentioned here that in addition to the baryon number conservation, the 
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parton densities (quark, anti-quark and gluon) also evolve satisfying the relation diri 1 = in 
absence of any particle production. Clearly, the evolution of the parton densities will also 
get modified due to the dissipative corrections. We estimate these dissipative corrections in 
a simple way from the following considerations. Since the shear viscosity scales as T 3 , while 
the bulk viscosity £ vanishes for a quark gluon plasma, one can define an acceptable range 
for the shear viscosity rj given by 0: 

2T 3 < V < Ut. (7) 

so that Navier-Stokes scaling theory is applicable to the expansion of the plasma. A baryon- 
free non-viscous plasma (/i=0, i]=0) leads to eT 4 ^ 3 =eoT^ 3 . Using the upper bound of vis- 
cosity as r)=er/4 and substituting e = 3P in Eq.(6), we get er=e T . This means that T 
scales as r~ 1//3 for i] = and as r~ 1//4 for i] approaching the upper bound. Therefore, it is 
reasonable to assume that T scales as r -7 for any value of r] in the range 1/3 < 7 < 1/4. 
Since e scales as r~ 47 , the LHS of Eq.(6) can be expressed in terms of 7 as, 

de e + P 4, N e 

7T + = -l-3 7 — 8 

or r 3 r 

On similar grounds, the evolution of parton density rii ( i=q, q, g ) can be described in one 
dimension as (since scales as t -37 ), 

^ + ^ = (l-3 7 )-. (9) 

OT T T 

For an ideal plasma, 7 = 1/3 and the RHS of Eqs. (8) and (9) vanish, while for other values 
of 7 the RHS is non-zero which can be interpreted as the correction due to dissipation. 
Comparing RHS of Eq.(8) with Eq.(6), we can identify 

77 = (1 — 37)er. 

It is also interesting to note that 77 = er/4 for 7 = 1/4 is consistent with the upper limit of 
viscosity. In terms of 77, Eq. (9) can be written as 

drij 7M _ ryrn 
dr r er 2 ' 

The dissipative correction appearing in Eq. (10) is derived for a baryon-free plasma. How- 
ever, as shown in Appendix A, the above correction is also valid for fi ^ if e and n q (q) are 
defined appropriately for a baryon-rich plasma. Finally, the evolution of baryon density can 
be obtained from Eqs. (9) and (10) as 

dn B n B n B r]n B 

~^— + = (1-37) = T- I 11 ) 
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where tib = n q — nq. It may be noted here that although the baryon conductivity vanishes 
for fi = 0, the dissipative corrections for the quark and anti-quark densities are finite as seen 
from Eq.(10). It may be mentioned here that for a baryon-symmetric matter, the energy 
momentum equation given by Eq.(6) is sufficient to describe the dynamical evolution of 
the plasma. However, Eq.(10) has been derived for consistency and will be useful to study 
the density evolution in presence of particle production particularly when the chemical 
equilibration is not complete. 

The dynamical evolution of a baryon-rich plasma can be studied by solving Eqs. (6) and 
(11) numerically. For the plasma phase, we consider only the u and d quarks for which the 
energy and baryon densities are given by 
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2 mQ/ x 



n B = "-T\x + ^). (13) 
where x = fi/T. The shear viscosity rj depends on both T and n in a complicated way ||11|| . 



However, as shown in Appendix B , we can approximate rj rjoe/T which scales as T 3 in 
case of a baryon free plasma. It is found that r/o weakly depends on T and /z. Therefore, we 
will treat it as a constant parameter for the subsequent calculations. 

Figures. 1(a) and 1(b) show the time evolution of the ratio of energy density e/eo and 
total entropy (ar)/(cr r ) for a plasma at x (/i/T)=1.0, energy density e = 9.0 GeV/fm 3 
(corresponding to an initial temperature of T = 0.258 GeV ) with different values of viscosity 
coefficients t]q= 0.0, 0.5 and 1.0 . As seen from the Figures 1(a) and (b), with increase 
in viscosity, the energy density evolves faster than that corresponding to the isentropic 
expansion (770 = 0) and also extra entropy is produced when the medium is viscous. Figure 
1(b) also shows that most of the entropy is produced in the first few fm/c and the asymptotic 
value of ar is approached quickly. So these results are consistent with the findings of |7j 
except that the present work has been extended to finite baryon density. As mentioned 
before, the shear viscosity may have an upper bound as given by Eq. (7). However, in the 
above calculations, we only carry out a parametric study to see the effect of viscosity on the 
evolution of the energy density and entropy production although some of the r/o values may 
exceed the upper limit. 

Figures 2(a) and 2(b) show the variation of T/T and <jt/(ct t ) as a function of r for 
different baryon densities of Xo=0.0, 1.0 and 2.0 at fixed 770 = 1 • and eo=9.0 GeV/fm 3 . For a 
given temperature T, the viscosity of the plasma at finite baryon density is higher than the 
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baryon-free case (since e increases with /i). This is true even if the energy density is held 
fixed since the plasma will have lower initial temperature at finite baryon density (resulting 
in higher viscosity). Since the viscosity increases with x, the temperature T evolves slower 
and more entropy is produced as compared to the baryon free case. It is also interesting to 
note that whatever be the value of viscosity, both T and fi scale with the same powerlaw as 
t~ 7 (see the solid circles in Figure 2a for [ij \Xq which overlap with the curve for T/T .) This 
means the ratio of \x and T (i.e. x) remains independent of r, consistent with the assumption 
made in appendix A. This is true even in case of an ideal (rjo = 0) baryon-rich plasma where 
both /i and T scale as r~l/3 i.e. both obey Bjorken's scaling. However, we have found 
that x decreases with r if the hydrodynamical expansion proceeds along with the chemical 
equilibration ||12|| . The r dependence of the chemical potential /i and the temperature T will 
be different due to parton production during chemical equilibration. 



III. EFFECT OF VISCOSITY ON CHEMICAL EQUILIBRATION FOR BARYON 

FREE PLASMA 

The plasma in equilibration means there is no net parton production and the momentum 
distribution of the partons can be described either by Fermi-Dirac (FD) or by Bose-Einstein 
(BE) distribution function. The picture emerging from the numerical simulations of parton 
cascades supports the formation of a plasma at RHIC and LHC energies which is initially 
gluon-rich. Rapid gluon scattering at a high initial temperature will lead to substantial 
thermalization of the plasma in a very short time while it will be far off from chemical equi- 
librium. Therefore, process of chemical equilibration will proceed along with the dynamical 
evolution of the plasma. The process of chemical equilibration will accelerate the cooling of 
the plasma further and the plasma may not achieve equilibrium by the time temperature 
drops to the critical value 0-§J. However, all theses calculations assume an ideal plasma 
where the effects of dissipative processes are neglected. The chemical evolution of a viscous 
plasma will be different than its ideal counterpart due to extra entropy production and 
viscous heating of the plasma. 

The distribution functions for quarks (anti-quarks) and gluons appropriate for a plasma 
which has achieved kinetic equilibrium but is far off from chemical equilibrium are given by 

f (-) = Mfi • f = ^ (141 

where Aj's give the measure of the deviation of the distribution functions from the equilibrium 
values. The chemical equilibrium is said to be achieved when Aj — > 1 and the baryo-chemical 
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potential /x attains an equilibrium value fj, eq consistent with the baryon current conservation. 
It may be mentioned here that although the plasma is dissipative, we use the FD and BE 
distribution functions for quarks, anti-quarks and gluons. However, the temperature and the 
chemical potential will evolve as per viscous hydrodynamics and will differ from their ideal 
counter parts. As discussed in Appendix A, the present formalism of viscous hydrodynamics 
is not valid for a baryon rich plasma undergoing chemical equilibration. Since x depends 
on r, additional corrections are needed in the rate equations describing the evolution of 
the viscous plasma. Therefore, in the following, we have studied the process of chemical 
equilibration for a viscous plasma only at zero baryon density. 

We restrict our considerations to the following two reaction mechanisms for the equili- 
bration of the parton flavours gg ggg and gg ^ qq. The evolution of the parton densities 
are governed by the master equations 

fy(n£) = (R 2 ^3 - #3^2) - (R 9 ^ q - R q ^ 9 ), (15) 

d^) = d IM (n^) = (R g ^-R^ g ). (16) 

where R 2 ^z and ^3^2 denote the rates for the process gg — > ggg and its reverse, and R g ^ q 
and Rq^g are for the process gg — > qq and its reverse respectively. Note that the density 
four vector n % = nu % + u l now includes a dissipative correction z/. The rate equations are 
in general quite complicated to solve. For baryon- free plasma (/i=0), the authors in [|J, 
have used a simple factorization for the RHS of the above equations based on a classical 
approximation i.e. using the Boltzmann distribution function for quark, anti-quark and 
gluon and eliminating the Pauli blocking and Bose enhancement factors in the final states. 
However, it was shown by us recently |l!J that the same factorization can be used, even by 



including full quantum statistics, particularly when the plasma is highly unsaturated. Using 
the viscous hydrodynamics as derived in the previous section and with the factorization 
given in Ref M, the Eqs. (15) and (16) can be written as 



^ + ^ = n 9 Rs(l - X g ) - 2n g R 2 (l -^) + M (17) 
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d ^ + ^ = n 9 R 2 (l-^h) + ^, (18) 



8r r y ^ X 2 g ' er 2 



dn- q n- q X q Xq nq V 
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The above set of equations differ from the corresponding equations of ref. [4] with the 
last terms containing the viscosity factors. The density weighted cross sections R3 and R 2 
are given by [[| 

#3 « 2.1 a 2 s T ^2X g - XI ; R 2 « .2AN f a 2 s X g T In(^). (20) 

Ol S Ag 

Using the distribution functions given by (14) for quarks, anti-quarks and the gluons, the 
energy and number densities can be written as: 

e = T 4 [a 2 X g + b 2 (X q + Xq)}. (21) 

with a 2 = 87r 2 /15; b 2 = Nf(77r 2 /40) ; where Nf is the dynamical quark flavours. 
Similarly, the number densities for gluon, quark and antiquark are: 

16 

n g = XgCuT 3 - ai = — ({3), (22) 

n q = W 3 ; b x = JL C(3)iV>, (23) 

n q = XgbiT 3 . (24) 

The above set of Eqs. (17-19) along with Eq. (6) were solved using fourth order Runga- 
Kutta method, to obtain the time dependence of the temperature T, the quark and anti- 
quark fugacities X g , X q (= Xq). The initial conditions were taken from the HIJING calcu- 
lations both for RHIC and LHC energies g. We take T = 0.57 GeV, A 90 =0.09, X q0 =0.02 
at r =0.31 fm for RHIC energy and T=0.83 GeV, X g0 =0.14, A g0 =0.03 and r =0.23 fm for 
LHC energy. Since, energy density depends on fugacities and chemical potential, viscosity 
will also depend on them implicitly. 

Figures (3) and (4) show the variation of T, X g and X q with proper time r for rj =0.0 
and 0.2. As mentioned above, these calculations correspond to a non-equilibrated plasma 
at zero baryon density. In the absence of dissipation, i.e. ^o=0, the plasma cools faster 
than that predicted by the Bjorken scaling (see the curves with solid circles in Figure 3(a) 
and 4(a) ), since additional energy is consumed in in the process of chemical equilibration. 
A further rise in equilibration rate (i.e rise in fugacities) will result in faster cooling of the 
plasma. The reverse happens in presence of viscosity. As seen from the Figures (3) and 
(4), during chemical equilibration, the quark and gluon equilibration rates become faster 
when the medium becomes more viscous. This, however, does not result in a faster cooling 
of the plasma as expected. So in the case of a viscous plasma, the chemical equilibration 
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becomes faster and the temperature drops more slowly due to generation of heat. Since the 
presence of viscosity makes the equilibration faster and also increases the life time of the 
plasma phase, it would be interesting to see whether chemical equilibration can be achieved 
when the medium is viscous. While this requires the precise knowledge of viscosity , the 
short-dashed curve in Figures (3) and (4) shows the result that can achieved with the upper 
limit of viscosity 77 = (er)/4 ( i.e 7] = ^f). These calculations show that even for the 
limiting case, the temperature may still drop faster than the Bjorken scaling and the plasma 
may not achieve chemical equilibration by the time T drops to T c « 0.2 GeV, provided the 
initial conditions taken from HIJING calculations are prevalent at RHIC and LHC energies. 
The above results are obtained for the chemical evolution of a viscous plasma at zero 



baryon density. In a recent paper, we have shown that |12| in case of a non-viscous plasma 
at finite baryon density, the rate of chemical equilibration slows down particularly for gluons 
resulting in slower cooling of the plasma. Therefore, in the presence of baryon density, the 
cooling rate of the viscous plasma is expected to slow down further. However, it is difficult 
to conclude about the net gluon equilibration rate since it is suppressed at finite baryon 
density. 



IV. SUMMARY AND CONCLUSION 



We have studied the properties of a viscous quark gluon plasma at finite baryon density. A 
set of equations that describes the space time evolution of the energy-momentum and parton 
densities has been derived based on a simple scaling law as well as by solving Boltzmann 
approximation under boost invariant scenario. The total entropy is conserved in case of an 
ideal plasma and the temperature and chemical potential evolve as per the Bjorken's scaling. 
However, for a viscous plasma additional entropy is produced due to viscosity both at zero 
and finite baryon density. In fact, finite baryon density makes the plasma more viscous as 
energy density increases with \i. Even at a fixed energy density, viscosity rises due to the 
lower initial temperature. Therefore, plasma at finite baryon density produces more entropy 
as compared to the baryon free case. We have also studied the chemical equilibration of the 
plasma at zero baryon density using the rate equations appropriate for a viscous plasma. 
It is found that due to viscosity , the temperature evolves more slowly, while the fugacities 
evolve faster than the ideal plasma. The slow cooling of the plasma and the fast equilibration 
rates of gluon and quark (anti-quark) will affect the thermal photon and di-lepton yields 
which are the important probes to study the dynamical evolution of the quark gluon plasma. 
The thermal open charm production will also be enhanced particularly at high Pj- or at high 
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invariant mass regions. The additional entropy which is generated due to viscosity will also 
reflect on an increased multiplicity distributions of the final particles. Further studies in 
this direction will be important to understand the behavior of the dynamical evolution of 
the quark gluon plasma produced in relativistic heavy ion collisions. 

To conclude, larger transverse momentum associated with the hydrodynamic expansion 
would be reduced as collective flow velocities are dissipated into heat. The dissipative effects 
may dampen the fluctuations which otherwise are expected to serve as signatures of unusual 
phenomena of QGP phase transitions. Since most proposed observables of the plasma are 
sensitive to the full space-time history of the reaction, dissipative phenomena must be taken 
into account if quantitative predictions are to be made. In future, we plan to estimate the 
dilepton and photon yields from QGP phase with the inclusion of these dissipative effects. 
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APPENDIX A 

The dissipative corrections r lk and v % in Eqs. (4) and (5) can be evaluated using the 
Boltzmann equation p^d^fa = C a (f), which describes the evolution of the Wigner densities 
f a (x,p) of the parton of type a in terms of the collision integrals C a . At equilibrium, 
the collision integral C(f) vanishes and the distribution function f (p,fi,T) becomes Bose- 
Einestein or Fermi-Dirac type. However, in the presence of dissipation, the hydrodynamical 
solution assumes only local equilibrium where 

f H = fa(p^(x),ti(x),T(x)). (A.l) 

with the four velocity u, chemical potential fi and temperature T being functions of x^. 
However, p^d^ju ^ 0, while C(fu) = 0. So fn cannot be a solution of the Boltzmann 
equation and there has to be a correction to fu by an amount 5f that is first order in 
gradients of u, /x and T. Writing / = f H ± 5f, the Boltzmann equation can be written as 

P< I dJ H =p a uJ-t. (A.2) 

where r c is the relaxation time which depends on time . In terms of 5f, the dissipative 
corrections are given by JF| 
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-ik 



(A.3) 



(2tt)^° 

The Boltzmann equation (A. 2) in (1+1) dimension has the simple form: 

dfn ^ df H 5f 



(A.4) 



dt 



dz 



(A.5) 



where v z is the parton velocity. As in |T3|, we further assume that the central rapidity regime 
is initially approximately Lorentz invariant under longitudinal boost so that the conditions 
at point z at time t are the same as those at z = at proper time r = (t 2 — z 2 ) 1 ^ 2 . Hence 
the left-hand side of Eq. ( A.5 ) at z = becomes simply the time derivative at constant p z t 

m 



df H (p±,Pz,t) 



dt 



5f(p±,Pz,t) 



(A.6) 



Therefore, the integral equations A(3) and A(4) can be evaluated using Sf given by 

df H (p±,Pz,t) 



Sf 



Tr- 



Ot 



Pzt 



f t 9 



(A.7) 



Pzt 



where /o = (1 =F fo) f° r fermions and bosons. An ideal plasma in one dimension expands 
following the Bjorken's scaling where the temperature T decreases as t _1//3 . In presence of 
dissipation, there is an extra entropy production and the plasma cools slowly. However, we 
can still assume a power law dependence for T oc t~ 7 where 1/4 < 7 < 1/3. Thus, the time 
derivative of (p/T) can be evaluated keeping p z t constant 



5/ = r c /o/o^(cos 2 ^- 7 ). 



(A.8) 



where p z = p cos 9 and p 2 = p\ + p 2 . Here we have assumed a baryon-free plasma ( 
fi = 0). However, for finite baryon density, p in Eq. (A.7) needs to be replaced by (p ± p) 
for anti-quark or quark distribution function respectively. Consequently, there will be an 
extra contribution arising from the term dx/dt where x = p/T. In case of a non- viscous 
plasma, by solving the equations for the conservation of entropy (at=const) and of baryon 
number (nst=const), it can be shown that both p and T decrease as t -1 / 3 even though the 
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plasma is baryon-rich. This implies that the ratio x remains independent of t. Therefore, it 
is reasonable to assume that x will also remain independent of t even for a viscous plasma 
at finite baryon density. This would mean both \x and T should follow the same power 
law as t~ 7 . In fact, this assumption is consistent with the results that we find in section 
II (see solid circles for /x//xo i n figure 2a). Since x does not depend on t, there would not 
be any extra correction and Eq.(A.8) can be used for Sf with distribution function f (for 
the quark and anti-quark) defined appropriately for a plasma at finite baryon density. The 
above assumption is not valid for a plasma undergoing chemical equilibration. It is shown 
in a recent work that x decreases with t due to parton production during chemical 



equilibration, although the time dependence can be neglected for small values of x. Since 
x varies with t, Eq. (A. 8) is not valid in case of a chemically equilibrating plasma at finite 
baryon density except when x=0 or very small. 

Using 5f given by Eq. A. 8 and BE or FD distribution for /o, the dissipative corrections 
for baryon-rich plasma, can be evaluated from Eqs. (A. 3) and (A. 4). Consider the case 
i = k = 0, for which r 00 and u° can be written as 

r 00 = Y t J 7^3 /V (/iW + f*M + f-J-^D d P J M («^ - 7)«n*, (A.9) 

u ° = TtJ^fJ p3 {fqfqT! ~ hI ^ ] dp I de {cos2e ~ l)sind - (A ' 10) 

Replacing the 9 integration by (2/3)(l — 37) and using 

POO _ poo 

/ P n fofodp = nT p^fodp, 
Jo Jo 

Eqs. (A.9) and (A. 10) can be written as 

^ 00 = ff(l-37h (A.11) 



i/ , = j(l-3r f )(n q -n q ). (A.12) 

where e = e q + e 5 + e g is the total energy density and n q , rig are the quark and anti- 
quark densities. Similarly, it can be shown that r 0z and v z will vanish due to vanishing 9 
integration. Recall that the terms which need to be added to the hydrodynamic Eqs. (3) are 
UkdiT %k and diV % . Since, we consider the expansion in (1+1) dimension around the central 
rapidity regime at z=0, the terms which are of interest to us are d^T 00 and Oqu . Although, 
gluons equilibrate much faster than the quarks and anti-quarks, we will assume an effective 
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relaxation time for gluon, quark and anti-quark given by r c 1 oc Ta 2 s ln(a s 1 ) fl4|| . Neglecting 
time dependence of a s one can have 

doT 00 = -^(3 7 +l)(l-3 7 )^, (A.13) 

^° = -(2 7 + l)(l-3 7 )||. (A.14) 

Therefore, in the presence of dissipation , terms proportional to e/(Tt 2 ) and n^/(Tt 2 ) are 
to be added to the Eq. (3) describing the evolution of the energy momentum and parton 
densities. It is also interesting to note that for an ideal plasma 7 = 1/3 so that the dissipative 
correction vanishes. However, comparing Eq. A(13) with RHS of Eq. (6) which has been 
derived from Eq.(4) directly, we can write 

dor 00 = (A.15) 

where r\ = (3 7 + 1)(1 — 3 7 )e/T. 

Similarly, in terms of rj, we can have 

d v = -c^r- (A.16) 

Therefore, we can write the following evolution equations in terms of r\ 

de e + P _4r] 



(A.17) 



dn B n B Cn B r] 

~W + T~—- (A ' 18) 

where C = (27 + l)/(37 + 1). It is also possible to write the density evolution for quarks 
and anti-quarks separately as: 

dn q n q _ Cn q r] 

dt + t et 2 ' 1 ' 

dn- q n q _ Cn q 7] 

¥ T" { } 

The above equations have the same form as that of Eqs. (6), (10) and (11). The value of 
C lies between 0.83 and 0.86 depending on the value of 7 . Ideally, we expect C=l from 
the requirement that in the absence of baryon density all the Eqs. A(17)-A(20) should give 
same solution for temperature. This deviation (C < 1) arises due to various approximations 
that are used in deriving 5f. For example, we have used an effective relaxation time r c both 
for quarks and gluons although r q > r g . Similarly, we use T = At' 1 where A is taken as 
constant, although it may depend on t. However, the deviation is not too big and we will 
use C=l for consistency in the subsequent studies. It is also important to note that the 
form of viscosity rj oc e/T is also consistent with the results of [IT] which was obtained using 
finite theory QCD (see Appendix B). 
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APPENDIX B 



The shear viscosity co-efficient using the relativistic kinetic theory for a massless QGP 
under relaxation time approximation can be written as JTIlinjJ 

4 
15 

where = q, q, g) is the energy density of particle type i and Aj is the mean free path, which 



is the inverse interaction rate, \ = 1/IY Recently, Hou and Jiarong |TT[ have evaluated the 
interaction rates for quarks and gluons as for baryon rich plasma using finite temperature 
QCD 

T q = Tu, T 9 = -Tv. (B.2) 

with 

167T 

where D is the gluon damping factor given by 



v = ~^a 2 s ( - + ^x 2 ) - In a s + In 4 , + a s + D \. (B.3) 



D~ , N ' a - N % 2 ). (B.4) 

Finally, the viscosity co-efficient rj is given by 

* = * = ^- < R5) 

Therefore, the \i dependence of the viscosity arises through e q and v. We have studied v 
as a function of /x at various temperatures T using the following expression for the running 
coupling constant: 

a(T ' ^ = (33 -2^) In ( T 2 ^ 2 y.6 22 y ( B - 6 ) 

where A s =0.1 GeV and Nf=2. It is found that v does not depend on /i very strongly Even 
the variation of v with temperature is also not very significant. The viscosity primarily 
depends on e/T where e is the total energy density of the plasma. Therefore, we can write 

V = VoTf- (B.7) 
where r] is treated as a constant although it may depend on T and \i rather weakly. 
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FIGURES 
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FIG. 1. (a) The ratio of the energy density e/eo , and (b)the ratio 
of the total entropy {(tt) / {uqTq) as a function of proper time r. 
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FIG. 2. (a) The ratio of the temperature T/Tq and (b) the ratio 
of total entropy (<tt)/(<7oto) as a function of proper time r. The solid 
circles are for n/fio as discussed in the text. 
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FIG. 3. The evolution of (a) temperatute T , (b)gluon 
fugacity A 9 and quark (anti-quark) fugacities X q as a func- 
tion of proper time r with initial conditions To=0.57 
GeV, to=0.31 fm, A 9 o= 0.09 and X q0 = A 50 =0.02 for 
r/o=0.0 and 0.2. The short-dashed curves represent the up- 
per limit when rj = ^. The curve with solid circles in (a) is 
the Bjorken's limit. 
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FIG. 4. Same as Figure 3, but for initial conditions 
at LHC, i.e T =0.83 GeV, r =0.23 fm, X g0 = 0.14 and 
\o = Ago=0.03. 
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